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1. Gamma-ray Transparency: A Convolution Product 

Gamma-ray cosmology is a line of research that aims to understand the propagation of TeV 
photons on cosmological distances, independently from the emission processes at play in the 
sources [1]. The main effect experienced by gamma rays on their journey is the production of 
electron-positron pairs in interactions with the extragalactic background light (EBL). The EBL is 
the second most intense cosmological background after the cosmic microwave background (CMB) 
and it encompasses the entire UV to far-infrared emission from stars and galaxies since the epoch 
of reionization. The H.E.S.S. and Femii-LP^ collaborations discovered the long-sought imprint of 
the EBE on the spectra of bright sources called blazars [2, 3] by jointly modeling the smooth intrin¬ 
sic spectral curvature expected from radiation processes in the sources together with absorption by 
the EBE. In [4], we extend such analysis methods and investigate the archival gamma-ray spectra 
of blazars published over the past twenty years. 

The fraction of gamma rays, exp(— t), surviving interactions with EBL photons is a function of 
the gamma-ray optical depth x{Eq,z), where Eq is the gamma-ray energy measured in the lab frame 
and z is the redshfit of the source [5]. The optical depth is the product of the distance between the 
emitter and the observer, the density of target photons, and the pair-production cross section. It thus 
results from a triple integration over the distance element, accounting for the cosmological model, ^ 
over the energy of the target photons, and over the incidence angle, with EBL photons assumed 
to be isotropically distributed in the comoving frame. We have shown in [4] that the integration 
over the incidence angle can be reduced analytically without any loss of generality. Going further 
requires an approximation, namely that the energy and redshift dependence of the EBL photon 
density can be decoupled, as initially proposed by [6]. The optical depth can be then computed as 
a linear function of energy, multiplied by the convolution product of the EBL specific intensity at 
z = 0, v/y, and of a function that integrates the cross section over the line of sight, which we name 
the “EBL kernel,” 
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where Oj is the Thomson cross section, and is the rest energy of an electron. The EBL 
intensity is here a function of ln(eo/'WeC^) = ln(/jv/mgC^), with eo the EBL photon energy in the 
lab frame and V the associated frequency, while the EBL kernel is a function of ln(£'o£o/”^eC^)^ 
and the convolution integrates over the dependence on Eq to leave a function of ln{Eo/m^c^). 

The EBL kernel, whose full expression can be found in [4], is defined above fhe pair-producfion 
fhreshold, (1 -\-z)^EQeQ/mlc‘^ > 1. The kernel shows a broad full widfh af half maximum, of abouf 
5 in fhreshold unifs. The peak posifion and amplifude are increasing funclions of redshiff, buf do 
nol sfrongly depend on fhe paramefer /evoi> which describes fhe evolution of fhe EBL wifh redshiff, 
as shown in Eig. 1. In [4], we funed fhe evolution parameter, /evoi = 1-7, on fhe models of [7] 
and [8], resulting in absorption deviations fo fhe models smaller fhan 15% up fo a redshiff z ~ 0.8. 
This is on fhe order of fhe sysfemafic uncerfainfy on fhe flux quoted by ground-based gamma-ray 
insfrumenfs, jusfifying our approach. We assume a single evolufion paramefer for fhe enfire EBL 


*We assume here a concordance cosmological model with Hq = 70km s * Mpc SIm = 0.3 and = 0.7, but 
other cosmological models could in principle be taken into account within our formalism. 
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Figure 1: EBL kernel for a source located at z = 0.1 and various evolutions. 

spectrum between ~ 0.1 and ~ 100/fm, but a natural extension of our approach could consist in 
developing a wavelength-dependent parametrization, possibly extending the applicability of this 
method beyond z = 1. 

2. Extracting the EBL Signal 

2.1 Dataset and Deconvolution Method 

We have collected published gamma-ray spectra from Whipple, HEGRA, CAT, Tibet, ARGO- 
YBJ, TACTIC, HESS, MAGIC, and VERITAS, with observations spanning the period 1995-2013. 
Each spectrum is described with either a power-law, a log-parabola, or an exponential cut-off 
power-law model, whichever is preferred by the ht accounting for EBE absorption. Among the 
106 spectra from 38 sources that we study, 86 spectra originate from known-redshift BE Eacertae 
blazars, which we use for “gamma-ray cosmology,” i.e. reconstruction of the EBE, constraints on 
the Hubble constant, and limits on Eorentz-invariance violation. We also collected the spectral 
characteristics of the sources observed by Fertni-LP^ at lower gamma-ray energies, when contem¬ 
poraneous with ground-based observations. Assuming that both the low- and high-energy gamma 
rays originate from the same radiative component, which is the current understanding of BE Eac¬ 
ertae emission, the Eerm/-EAT observations below lOOGeV constrain the unabsorbed part of the 
spectrum, imposing a maximum hardness on the intrinsic emission at higher energies. 

Eor htting purposes, we bin the EBE specihc intensity in the wavelength range 0.26 — 105 /im 
in eight bins of width Ain A = 0.75. The maximum wavelength is determined by the pair-production 
threshold for the highest-energy gamma-ray spectral points, while the number of bins and minimum 
wavelength are limited by the gamma-ray statistics. The overall EBE spectrum is represented by a 
sum of Gaussians centered on each bin, with hxed widths and positions, and with free normaliza¬ 
tions that are the ht parameters. The parameters of each intrinsic gamma-ray spectrum are left free 
to vary during the minimization procedure over the EBE parameters, following the approach of [2]. 
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In addition to the quality of the fit of the gamma-ray spectra and to the maximum hardness 
limits, we can further constrain the EBL parameters using results from direct measurements and 
galaxy counts, which we largely extracted from [9]. The former are sometimes prone to contam¬ 
ination from local foregrounds, such as zodiacal light, and are considered as upper limits in our 
minimization procedure. The latter only account for resolved populations of sources, and not for 
any truly diffuse component or undetected classes of objects such as primordial stars. We thus 
consider the results from galaxy counts as lower limits on the EBL intensity. More details on the 
data, deconvolution method, and minimization procedure can be found in [4]. Eor those gamma- 
ray spectra where the original authors have agreed, we plan on submitting the data points we have 
assembled to the ASDC SED builder and TeVCat,^ so they will be easily available to others. 

2.2 Best-fit EBL Spectrum 

The eight-point best-fit spectrum resulting from the 86 gamma-ray spectra and from local 
measurements is shown in Eig. 2. Eixing the choice of intrinsic spectral model of each spectrum 
and leaving their parameters free, the best-fit EBL spectrum not accounting for local measurements 
is preferred by a likelihood ratio test at the 11a level with respect to no EBL. Eor comparison, 
we also show in Eig. 2 the result of a model-dependent approach where the normalization of the 
EBL spectrum modeled by [8] is left free to vary, reaching here a best-fit value of 1.13 ±0.07. 
We find that state-of-the-art EBL models [7, 8, 10] represent the data equally well as our model- 
independent approach based on Gaussian sums, and we estimate the total brightness of the EBL in 
0.1 — 1000 Jim as 62± 12nW m^^ sr^', or 6.5 ± 1.2% of the brightness of the CMB. 

Mild discrepancies between the model-dependant and -independent approaches are observed 
around 10/im, which could possibly trace the signature of polycyclic aromatic hydrocarbons. We 
are in the process of investigating the impact of different reconstruction techniques, such as spline 
functions, to assess the robustness of these features. Erom the difference between gamma-ray 
constraints and galaxy counts, we can exclude any contribution to the EBL intensity larger than 
~ 6nW m^^ sr^' from unresolved sources between 0.1 and 10/im, which rules out reionization 
scenarios involving miniquasars or primordial stars such as presented in [11]. 

3. Beyond the Standard Models of Cosmology and Particle Physics 

3.1 Hubble Constant 

The Hubble constant traces the expansion rate of the Universe. Measurements based on the 
cosmic distance ladder point to an expansion rate Hq = 72.5 ± 2.5km s^' Mpc^^, in ~ 2a ten¬ 
sion with the latest CMB constraint from Planck, Hq = 67.0 ± 1.2km s^^ Mpc^^ [12]. Though 
this mild discrepancy might be merely of statistical origin, further tensions between local and dis¬ 
tant measurements of the Hubble constant might reveal some limitations in the standard model of 
cosmology. 

A number of authors, e.g. [13, 14, 15], proposed to use gamma-ray absorption as a probe of 
the Hubble constant. The first-order dependence of the optical depth on Hq is shown in Eq. 1.1 and 
was used directly by [14] to set a lower limit on the Hubble constant. Similarly, by studying the 

^http://tools.asdc.asi.it/, http://tevcat.uchicago.edu/ 
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Figure 2: Specific intensity of the EBL. Galaxy counts and direct measurements are extracted from the liter¬ 
ature and represented by upward- and downward-going arrows, respectively. The model-independent result 
from this study based on gamma-ray observations is shown with blue points, while the model-dependent 
result based the EBL intensity of [8] is shown with a filled gray region (la confidence level). 

ratio of the EBL reconstructed from gamma-ray observation and that inferred from galaxy counts, 
we reconstruct //q = 88 i 8stat i 13sys kni s^' Mpc^\ for a fixed evolution parameter /evoi = l-V. 
Due to its large uncertainties, our measurement is compatible at the 1.4 a level with other state- 
of-the-art constraints, so that no hint for physics beyond the standard model of cosmology can be 
claimed at this stage. The authors of [15] notice that the density of photons could also depend on Hq 
for a given set of galaxy observations. Within our formalism, this would translate into varying the 
evolution parameter as a function of the Hubble constant. Marginalizing over 0.5 < /evoi < 2.5, we 
reconstruct instead //q = 88 ± 13stat =t 13sys km s^' Mpc^\ which does not affect our conclusions. 

3.2 Lorentz-Invariance Violation 

The unification of quantum field theory, which describes physics at the particle scale, and of 
general relativity, relevant at cosmological scales or in extreme astrophysical conditions, is one of 
the most important quests of modern-day science. The search for a quantum-gravity theory has led 
to the development of new mathematical abstractions, such as loop quantum gravity or string theory 
(see, e.g., discussions in [16]), which have nonetheless sometimes been criticized for their lack of 
testable predictions. Phenomenologists have to a large extent deflected the criticism and bridged 
the gap between theory and experimental physics, by studying generic leading-order corrections to 
“classical” theories, and their testability in fields such as cosmology, particle physics, astroparticle 
physics, and gamma-ray astronomy [17]. 

One of these effects originates from a UV correction to the dispersion relation of particles, 
yielding a norm of the momentum four-vector that is not a Lorentz invariant [18]. Assuming that the 
momentum four-vector is conserved between the initial and final states, one can compute a modified 
pair-production thr'eshold. In the case of a subluminal first-order correction, the threshold condition 
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Figure 3; Gamma-ray transparency as a function of energy in the classical case (solid lines) and with a 
modification of the pair-production threshold at the Planck scale (dashed lines). 

becomes (1 -|-z)^£'o£o/m^c‘^ > 1 + {Eq/Ejus)^, with Ejun = ~ 29.4TeV x 

(£'QG/£'pianck)^'^^ and where Eqq is the quantum gravity energy scale, assumed here to be on the 
order of the Planck energy scale Spianck = y/hcfjG = 1.22 x 10^^ eV. Interestingly, despite a 
difference of fifteen orders of magnitude between gamma rays observed on Earth and the Planck 
energy scale, the pair-production threshold is already affected at the ~ 10% level around 15TeV. 

Following the work of [19], we account for the modified threshold with a formalism similar 
to that presented in Eq. 1.1. The net effect is an increased transparency to gamma-rays above 
lOTeV, largely independent of redshift, as shown in Fig. 3. We constrain the effect by fitting 
to the 86 spectra a model including the quantum-gravity energy scale as a free parameter. The 
EBE parameters and the intrinsic parameters are left free during the minimization procedure. The 
likelihood profile as a function of Eqq shows a maximum around the Planck energy scale, though 
only significant at the 2.4 a level. We set the first quantitative limit on the effect, accounting for the 
uncertainties on the EBE and on the intrinsic emission of BE Eac objects, as Eqq > 0.6 x Epianck at 
the 95 % confidence level. Future observations with instruments sensitive to higher-energy gamma 
rays, such as CTA [20], HAWC [21], and EHAASO [22], will be able to rule out or unveil a 
modification of the pair-production threshold at the Planck scale. 

3.3 Signatures of Ultra-High Energy Cosmic Rays and Axion-Like Particles 

The study of Eorentz-invariance violation is a natural, one-parameter extension of the standard 
picture, in which gamma rays are emitted at the source and pair produce on the photons emitted by 
stars and galaxies since the end of the dark ages. More complex extensions of the standard picture 
have been proposed. A first one, based on classical physics only, relies on the acceleration and 
escape of a significant number of ultra-high energy hadrons from the sources. These hadrons would 
then produce neutral pions in interactions with the EBE, generating the VHE radiation observed 
on Earth [23]. A second one involves the conversion of gamma rays to hypothetical axion-like 
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particles in magnetic fields of the blazar jet, galaxy, or galaxy cluster (e.g. [24]). These axion-like 
particles then propagate freely through the EBL and potentially convert back to gamma rays in the 
magnetic field of the Milky Way. Both effects have been claimed to explain some features of the 
experimental data, such as possibly harder-than-expected spectra at the highest energies or optical 
depths (e.g. [25]). 

The parameters of these models vary from source to source (e.g. strength of the magnetic 
fields, fraction of hadrons escaping the source, etc), so that detailed constraints require studies 
beyond the work presented in [4]. We can nonetheless question the need for such scenarios by 
investigating the match between the observed gamma-ray spectra and EBL absorption alone. Note 
that for sources with underconstrained redshifts, we fix the distances to their best estimates based 
on the various constraints found in the literature. Based on the results presented in Sec. 2.2, we 
reach the following conclusions: I) The 106 spectra we studied are well matched by gamma-ray 
absorption, with probabilities larger than 9 % for all of them. II) We do not find any tension with 
the hardness limit imposed by the Fermi-LKV observations. Ill) The level of EBL reconstructed 
from gamma-ray observations only is in good agreement with constraints from direct measure¬ 
ments and galaxy counts. IV) We do not find any correlation between the fit residuals and the 
energy or the optical depth, at odds with the conclusions of [25]. We conclude that the standard 
picture is, at this stage, sufficient to explain gamma-ray observations of blazars, without the need 
for scenarios involving axion-like particles or ultra-high energy cosmic rays. We encourage further 
constraints on these scenarios using state-of-the-art EBL models and a large unbiased sample of 
gamma-ray spectra. 

We gratefully acknowledge support from the U.S. National Science Eoundation award PHY- 
1229792. We would also like to thank Matthieu Bethermin, Michael Dine, Alberto Dominguez, 
Steve Eegan, Amy Eurniss, Joel Primack, and Stefano Profumo for helpful discussions. 
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